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ABSTRACT 
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of  an  existing  continuous  control  system  has  been  expanded  by  Professor 
Benjamin  C.  Kuo  of  the  University  of  Illinois  to  enable  digital  redesign 
of  a  forced  control  system.  This  report  presents  a  brief  description  of 
the  newly  expanded  technique  and  includes  an  example  of  its  application. 
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1.  Introduction 


A  technique  under  development  by  Professor  Benjamin  C.  Kuo 
of  the  University  of  Illinois  has  been  described  and  applied  tc  several 

design  examples.***'  The  technique  enables  the  designer  to  perform  a  digi¬ 
tal  redesign  of  an  existing  continuous  control  system.  The  technique 
was  originally  developed  for  application  to  certain  linear  regulator 
systems.  Recently  Professor  Kuo  expanded  the  capability  of  the  tech¬ 
nique  to  enable  digital  redesign  of  a  forced  control  system  which  enables 
the  technique  to  be  more  appropriately  applied  to  the  models.  This 
report  describes  briefly  the  extension  of  the  technique  introduced  and 
includes  an  example  of  its  application. 


2.  Theoretical  Basis 

Assume  a  linear  continuous  control  system  is  df scribed  l?  the 
state  variable  equations, 


x(t)  =  Ax(t)  +  Bm(t)  (1) 

and 

m(t)  =  'il(O)  r(t)  -  G(0)  x(t)  (2) 

where  x(t)  is  a  vector  representing  the  state  of  the  system;  A,  B,  G(0, , 
and  E(0)  represent  time- invariant  matrices;  andr(t)  is  a  vector  repre¬ 
senting  the  input  to  the  system. 

The  system  may  be  represented  by  the  block  diagram  shown  in  Figure  1. 
The  objective  of  this  investigation  is  to  find  a  saripled-data  system 
such  that  its  response  co  an  arbitrary  input  will  be  identical  to  the 
response  of  the  continuous  control  system  at  the  sum;: ling  instnts. 

Suppose  it  is  desired  to  replace  the  continuous  control  system 
[represented  by  Figure  1  and  Equations  (1)  and  (2)]  with  a  digital  ver¬ 
sion.  This  version  is  achieved  by  adding  a  sampling  device  ^vith  sample 
periol  T)  followed  by  a  zero  order  hold  (also  called  a  clamp  or  boxca\ 
device).  Such  a  control  system  is  represented  schematically  in  the  slock 
diagram  of  Figure  2  where  m(kT)  represeucs  the  output  of  the  zero  order 
hold  (shown  as  the  block  labeled  ZOH)  and  k  is  an  integer. 


Presented  by  Professor  Kuo  at  a  seminar  in  the  Astrionics  Labora¬ 
tory,  Marshall  Space  Flight  Center,  MASA,  Huntsville,  Alabama,  20-24 
July  1970  (modified  in  subsequent  correspondence  between  Professor  Kuo 
and  Dr.  Seltzer). 

2Seltzer,  S.  M.,  Digital  Redesign  of  a  Rocket  Attitude  Control  System. 
U.S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama,  January  1971, 
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The  goal  of  this  approach  is  to  make  the  sampled-data  system  of 
Figure  2  equivalent  to  the  continuous  system  of  Figure  1  in  the  sense 
that  m(t)  will  be  equal  to  m(kT),  where  t  lies  between  kT  and  (k  +  1)T, 
and  the  state  trajectories  of  the  two  systems,  represented  by  the  state 
Equations  (1)  and  (3),  are  similar  (in  the  sense  implied  above).  Then 
the  discrete  data  system  may  be  represented  analytically  by 


x(t)  =  Ax(t)  +  Bta(t) 


and 


m(kT)  =  E(T)  r(t)  -  G(T)  x(kT) 


(3) 

(4) 


The  state  vector  x(t)  is  used  rather  than  x(t)  to  denote  that  the  formerly 
continuous  system  now  is  sampled  and  has  a  different  feedback  matrix, 

G(T),  and  input  matrix,  E(T). 


The  development  presented  under  Section  2  of  Reference  2  remains 
initially  unchanged,  i.e.,  the  following  relationships  still  i  old  ior 
G(T): 


G(T)  =  G(0)  +  T 


ST 


T:.0 


(5) 


where 


=  \  G(0) [A  -  B  G(0) ] 

0 


(6) 


An  expression  for  E(T")  is  found  in  a  similar  manner,  leading  to  the 
equation, 


E(T)  =  E(0)  +  T 


BE  (T) 

St 


Jt  =  o 


(7) 


whe  re 


=  -  \  G(0)  BE (0) 


SE(T) 


(8) 


If  G(0),  E(0),  A,  and  B  are  specified  for  the  continuous  model,  G(T)  may 
be  obtained  for  the  sampled-data  model  by  applying  Equations  (5)  through 
(8). 


3.  Example 

A  primitive  linearized  model  of  a  rocket  dynamical  continuous 
control  system  may  be  characterized  by  the  state  Equations  (1)  and  (2) 
where 


•  T 

x  =  [0,  0] 


(9) 


A  = 


"0 

-0 


(10) 


B  =  (0,  c]T  ,  (11) 

E(0)  =  1  ,  (12) 

and 

r(t)  =  X,  t  2  0 

=  0,  t  <  0  .  (13) 

The  corresponding  block  diagram  representing  this  simple  time -invariant 
system  is  shown  in  Figure  3.  For  a  particular  analytical  model  of  a 
rocket  system  being  studied  by  the  U.S.  Army  Missile  Command,  a  value  of 
c  corresponding  to  passage  of  the  vehicle  through  maximum  dynamic  pres¬ 
sure  conditions  would  be  90  (in  units  consistent  with  the  other  terms 
used  :n  Equation  (1)].  To  obtain  values  for  control  gains  and  a^ 

that  will  yield  desirable  response  characteristics,  the  characteristic 
equation  corresponding  to  this  system  may  be  obtained  from  the  relation, 


where 


(14) 


x  =  A  x  +  B  X 


(15) 
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and 


(16) 


Equation  (14)  leads  to  the  characteristic  equation, 

2 

s  +  a1c  s  +  aQc  =  0  ,  (17) 

where  cha  natural  frequency  w  and  damping  ratio  £  may  be  written  as 


and 


(18) 


(19) 


If  a  value  of  6  radians/second  is  chosen  for  m  an!  1/./2  for  f 

n  ' v  1 


0(0)  =  [aQ,  a1'|  =  f 0.398,  0.094] 


(20) 


Using  these  values  for  G(0),  the  expected  response  to  a  step  input  X  may 
be  written  as 


ml  „  x-1 


sis  +  aLc  s  +  aQcl 


=  3.54  e‘4,25tcos(ut+ 135  deg) +  2.5  (21; 

and  is  plotted  in  Figure  4.  Application  of  the  final  value  theorem  leads 
co  the  steady-state  value  of  X/aQ  for  9(t) , 
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Applying  the  technique  described  in  Paragraph  2,  the  continuous 
system  of  Figure  3  may  be  converted  to  a  sampled-data  system  as  indicated 
in  Figure  5.  Equations  (3)  and  (4)  may  be  used  to  describe  the  system 
and  Equations  (6),  (8),  (10),  (11),  (12),  and  (20)  used  to  obtain  the 
expressions, 


dG(T)' 

3T 


f  G(0)  [A  -  BG(0) ] 


-  I  (-ao  ai)  •  (t  ■  ai) 


-  - 1  G<°)  BE(°)  -  -  i  v 

_T  =  0 


Using  Equations  (5)  and  (22),  ar.  expression  for  the  feedback  gains  may 
be  written: 


G(T)  =  a0«  al  “  aQ 


Similarly,  using  Equations  ('/)  and  (23),  an  expression  tor  E(T)  may  be 
obtained:. 


E(T)  =  1  -  |  a^T  . 


Using  Figure  5  and  Equations  (24)  and  (25),  the  system  transfer  function 
may  be  found : 


ilzl  _ 

X(z) 


“2  E(T) 


1  +  Z 


:)(7)(w) 


c  T  E  (T)  (z  +  1) 


2,2  +(a„c  T2  +  2\c  T  -  4),  +(a0c  I2  -  2^0  T  +  l) 


ESWTCT*'  W***  { ■--' m-  s*> V 


Application  of  the  final  value  theorem  to  Equation  (26)  leads  to  the 
steady-state  value  of  X/a_  for  9  when  the  system  input  is  a  step  of 
magnitude  X. 

If  the  numerical  values  of  Equation  (20)  are  used  for  a~  and  a  and 
c  is  set  at  90,  G(T)  and  E(T)  are  found  to  be: 

G(T)  «  [(0.398  -  1.692  T),  (0.094  -  0.2  T)]  =  ^aQ,  j  ,  (27) 

and 

E(T)  -  1  -  4.2488  T  .  (28) 


To  preclude  positive  feedback,  T  must  have  a  numerical  value  less  than 
0.398/1.692  =  0.24. 

Using  the  results  of  Equations  (27)  and  (28),  Equation  (26)  may  be 
rewritten  as 


aw  w*1  («(■))  -z'1 


2«2  +  [(39. a  -  152.3T)T2  +  (16.9  -  36.01)1  -  *]«  ♦  [(39. »  -  152.3r)I2  -  (16.)  -  36.01)1  +  2] 


(29) 


For  a  step  input  of  magnitude  X, 


X(z) 


_ Xz _ 

(z  -  1) 


(30) 


Values  of  0(nT)  may  be  found  from  Equations  (29)  and  (30)  for  various 
selected  sampling  periods,  T  (Figure  4). 


4.  Conclusions 

The  technique  for  perfjrming  digital  redesign  of  a  missile 
attitude  control  system  that  can  be  formulated  in  the  form  of  Equations 

(1)  and  (2)  is  simple  and  easy  to  apply.  As  the  chosen  value  for  T 


b 


approaches  the  limiting  value  of  0.24,  the  step  input  response  (0) 
becujieo  taui.e  sluggish,  becoming  unstable  when  T  >  0.24  (Figure  l). 


Figure  \.  Continuous  Model 


Figure  3.  Simplified  Rocket  Continuous  Cou'.rol  System 
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